Structural phase transitions in Ruddlesden-Popper phases of strontium titanate: 
ab initio and inhomogeneous Ginzburg-Landau approaches 
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We present the first systematic ab initio study of anti-ferrodistortive (AFD) order in Ruddlesden- 
Popper (RP) phases of strontium titanate, Sri+ n Ti n 03 n +i, as a function of both compressive epitax- 
ial strain and phase number n. We find all RP phases to exhibit AFD order under a significant range 
of strains, recovering the bulk AFD order as ~ 1/n 2 . A Ginzburg-Landau Hamiltonian generalized 
to include inter-octahedral interactions reproduces our ab initio results well, opening a pathway to 
understanding other nanostructured perovskite systems. 

PACS numbers: 73.21.Cd, 68.35.bg, 68.35.Rh 



Superlattices originating from various oxide per- 
ovskites are of great interest due to their rich properties. 
As examples, S^RuC^ of the n = 1 Ruddlesden-Popper 
(RP) family exhibits unconventional superconductivity 1 , 
ferroelectricity in multicomponent superlattices made of 
different perovskites (BaTi0 3 , SrTi0 3 and/or CaTi0 3 ) 
can be tuned by controlling the mixing ratio 2,3 , and, 
recently, even a superconducting two-dimensional elec- 
tron gas has been observed at the interface of the 
SrTi0 3 /LaA10 3 superlatticel 

The Ruddlesden-Popper series of structures, 
A n+1 B n 3n+1 = (AB0 3 )„(AO) (or, "RP-n") for 
n = 1,2,..., is the simplest prototype of such nano- 
ordered superlattices, consisting of a set of AO stacking 
fault planes interspersed regularly between every n 
layers of bulk perovskite octahedral. The RP-n se- 
ries most successfully synthesized to date is that of 
strontium titanate, with phases up to n = 5 reported 
to have been grown controllably with molecular beam 
epitaxy(MBE) 6 . The n = oo end member of the series, 
bulk SrTi0 3 , is also of great interest in its own right, 
with highly tunable properties such as strain-switchable 
ferroelectricity^, and even metallic and superconducting 
behavior under doping with oxygen vacancies^. 

Despite its technological and scientific potential, the 
RP-n series of strontium titanate remains relatively un- 
explored. The bulk material exhibits a rich strain- 
temperature structural phase diagram including both 
anti-ferrodistortive (AFD) and ferroelectric (FE) order- 
ing and various combinations thereof!- , with non-FE 
AFD order present at zero strain and temperature. In 
contrast, previous ab initio work has found there to be 
no AFD order at zero strain and temperature in the RP- 
n series for n — 1,2, with the expected bulk-like order 
only occurring in n = 3, 4, 5 n i 12 . However, the associ- 
ated phase transitions with strain, the underlying driving 
forces, the scaling of the behavior with n, and the nature 
of the large-n limit remains unexplored. 

This work addresses the above open issues by employ- 
ing the ab initio approach to study the AFD rotational 
transition as a function of both distance between stacking 
faults n and epitaxial strain (biaxial strain applied in the 
plane of the SrO stacking faults) . Such epitaxial strain is 



now accessible experimentally, with substrates available 
to provide a substantial range of strains to tune the ma- 
terial properties of strontium-titanate structure o 10 i 13 ' 14 . 
For this initial study, we here focus on compressive 
strains, which favor AFD order with the rotation axis 
perpendicular to the stacking fault planes. This allows 
us to investigate the transition to bulk behavior with in- 
creasing n, while freeing us from having to consider the 
more complex configurations associated with the multi- 
ple, in-plane rotational axes favored by tensile strains. 

Computational details - All ab initio calculations be- 
low employ the density-functional theory framework 15 
within the local-density approximation (LDA) and rep- 
resent the ionic cores with norm-conserving Kleinman- 
Bylander pseudopotential$i£. A plane-wave basis with a 
cutoff energy of 30 H expands the Kohn-Sham orbitals, 
and Monkhorst-Pack meshes^ of 2x2x2 or 2x2x1 (for 
RP-n structures of n < 2 or n > 3, respectively) sample 
the Brillouin zone. Finally, the resulting energy func- 
tional is minimized by the analytically continued conju- 
gate gradient method^ within the DFT++ software^. 

We employ supercells composed of two bulk regions, 
containing n strontium titanate bulk layers each, along 
with two extra SrO stacking fault layers. To properly 
represent the AFD order, each bulk layer contains two 
structural units of strontium titanate. To represent epi- 
taxial strain, we fix the in-plane lattice parameter a and 
relax the out-of-plane lattice constant c. All calcula- 
tions minimize the electronic wave functions to within 
10 /iH of the Born-Oppenheimer surface, the ionic coordi- 
nates until all forces on individual atoms |F| are less than 
0.1 mH/B, and the out-of-plane strain until it is deter- 
mined to within ±0.2% (1 H«27.21 eV, 1 B«0.0529 nm). 

With the ultimate goal of comparing ab initio re- 
sults to experiments on epitaxially strained thin films, 
we define epitaxial strains as e = (a — clb)/<ib, where 
as = 3.8376 A is the equilibrium lattice constant for 
cubic (non-rotated) bulk SrTi0 3 within our computa- 
tional framework. We find that small-n phases under 
sufficiently low compressive strains either exhibit AFD 
order with an out-of-plane rotation axis or remain non- 
ordered (centrosymmetric) . To obtain reliable transition 
points for low-n phases and to extract information about 
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FIG. 1: (Color online) Structural phase diagram of 
Sr n +iTi n 03„+i as a function of epitaxial strain and n: in- 
plane FE AFD order (x), out-of-plane FE AFD order (o), 
out-of-plane non-FE AFD order (+), non-ordered (*), non- 
FE AFD to non-ordered phase boundary studied in the text 
(solid line), additional phase boundaries (dashed lines). 

bulk SrTiC>3 that can be connected directly to our re- 
sults at low n, we sometimes enter regions of phase space 
where the material exhibits other instabilities. Figure Q] 
summarizes the phase space we considered and the phases 
we found in each region. Under sufficiently high compres- 
sive or tensile strains, AFD and FE orders coexist in the 
large-n phases with the order parameters oriented along 
either the out-of- or in- plane directions, respectively, just 
as in bulk strontium titanatei^. 

Concerned here only with the transition between the 
out-of-plane non-FE AFD phase and the non-ordered 
phase, for the calculations below, we suppress FE order 
and in-plane AFD order by imposing both C2 and mir- 
ror symmetry about the out-of-plane direction. Figure [1] 
shows that, for n < 4 and low strains, the ground-state 
phases are free of these suppressed orders, so these sym- 
metry restriction in no way affect our conclusions for the 
n < 4 phases. For the n > 5 phases, the in-plane rotated 
order will likely become more thermodynamically stable 
near the phase-transition point of interest to us. Thus, 
for large n, the rotated phases which we study may be 
only metastable near the transition point. Because we 
here are concerned with the nature of the non-FE AFD 
to non-ordered phase boundary and how the spacing n 
between stacking faults affect this transition, the impo- 
sition of the above symmetry restrictions again does not 
affect the central conclusions below. 

SV^Tifti, n = 1 phase - As a first look at the AFD 
phase transition under strain, we begin by noting that, by 
the symmetries of the 00c~ AFD ordering in the Glazer 
notation- , a single rotation angle 9 serves to define the 
state of rotation for each layer of octahedra. In a very 
simple picture in which the rotations in all layers in each 
bulk region are the same (as is the case for the n = 1, 2 
phases), we would have a single order parameter and the 
usual Ginzburg-Landau form for the free energy, 

F(9) = ^A9 i + a(e~e c )9 2 + F , (1) 
where both A and a are positive, and e c is the criti- 



FIG. 2: (Color online) Square of rotational order parameter, 
6 2 , versus epitaxial strain e (upper panel); square-root of tran- 
sition energy, V AE, versus epitaxial strain e (lower panel) : ab 
initio results (crosses), linear fits (dashed lines). 



cal strain of the transition, with the well-known solution 
that the ordered phase obtains for e < e c , with a nonzero 
rotation angle 9 2 {e) = a(e c — e)/A and an energy dif- 
ference relative to the centrosymmetric state (9 = 0) of 
^JAE{e) = a(e c — e)/\/2A. Thus, the square of rotation 
angle and the square-root of the energy difference should 
both vanish linearly as e approaches e c from below. (Here 
and below we regard compressive strains as e < 0.) 

Figure [2] shows that our full ab initio results for the 
n = 1 RP phase exhibit quite clearly these classic signa- 
tures of a second-order phase transition. The extrapo- 
lated values for the critical strain from the two observ- 
ables, 9 2 and V AE, are in surprisingly good agreement 
for such a simple model, -2.2% and -2.8%, respectively. 
In either case, it is clear that the critical strain is sig- 
nificantly below zero, consistent with the lack of obser- 
vation of AFD rotational order in the n — 1 phase near 
zero strain in previous ab initio workii. As we now show, 
the critical strain passes through zero and becomes pos- 
itive with increasing n, eventually approaching a value 
associated with the bulk material. 

iSr Jl+ i Tin Osn+i, n > 1 phase - For n > 1, multiple 
perovskite layers exist between faults and different ro- 
tation angles become possible for each layer (except for 
the n = 2 phase where the two layers are identical by 
symmetry). The upper panel of Figure [3] displays the 
rotation angles in each layer for the n = 5 phase at zero 
strain. These data exhibit the general trend in all of 
our data, namely that the rotation angles are smaller 
near the SrO stacking faults and approach the expected 
bulk value toward the center of the bulk regions. Associ- 
ated with the reduced rotation angle at the interface, we 
find an enhanced stretching of the octahedra neighboring 
the stacking faults resulting in rumpling of the SrO lay- 
ers. We also find off-center motion of the titanium ions 
of these boundary octahedra toward the fault layers, so 
that each region of bulk material, and thus the entire 
structure, exhibits no net ionic electric dipole moment. 
This suggests that the structure near the SrO stacking 
faults is such that the boundary octahedra prefer stretch- 
ing and formation of a ionic electric dipole moment over 
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FIG. 3: (Color online) Ab initio (+) and Ginzburg-Landau 
model (o) results for layer-dependent rotation angles 6i versus 
layer position for the n — 5 phase at e = 0.0% (upper panel), 
and critical strain e c versus inverse-square of RP phase num- 
ber n (lower panel). Lines in upper panel are guides to the 
eyes; solid line in lower panel is best fit to ab initio data. 

rotation as a way to relieve compressive stress, ultimately 
suppressing the rotational instability near the faults. 

To locate the AFD transition, for the order parameter 
9, we take the rotation angle of the central one or two 
(if n is odd or even, respectively) perovskite layers, and 
extrapolate the linear behavior of 8 2 (which shows less 
numerical scatter than that for V AE), identifying the 
horizontal intercept as the critical strain 2 ' 1 -. The lower 
panel of Figure [3] displays the resulting critical strains 
for n > 2 as a function of 1/n 2 . With increasing n, the 
critical strain steadily increases toward the bulk value as 
the RP-n phases become more like bulk. The critical 
strains are negative for n = 1,2 and positive for n > 
3, consistent with and explaining the previous ab initio 
observation a 1 L 1 2 that only n > 3 phases exhibit AFD 
order when grown on a bulk SrTiC>3 substrate(e = 0%). 
Finally, the clear linear behavior in the plot indicates that 
the critical strain approaches the bulk value as 1/n 2 . 

To explore the origin of the above behaviors, we gener- 
alize the above Ginzburg-Landau free-energy model to in- 
homogeneous systems by including both inter-octahedral 
interactions and a spatially-dependent tendency toward 
rotational instability, 

nm) = E {\ A6 i + a ( e - + W* + *-i) 2 } . 

(2) 

where 0j is the octahedral rotation angle of the i-th layer, 
A and a are defined as above (Eq. [1} , e, represents an 
effective critical strain for the i-th layer, and the third, 



inter-layer term imposes the alternating sign ordering of 
0i associated with the Glazer 00c~ order of the material. 
To reflect the lessened tendency towards rotation at the 
interfaces, ej should have more negative values near the 
stacking faults, eventually approaching the critical value 
for the bulk material e c (n = oo) in the center of the bulk 
regions. In this initial work, we take £j to be the bulk 
value e c (oo) for all layers i except those neighboring the 
stacking faults, for which we will take a different value, 
6j = e, to be determined below. The standard stability 
analysis for this free energy is that there will be a second- 
order phase transition when the Hessian, evaluated at 
9i = for all i, ceases to be positive definite. 

To begin our analysis of the above free-energy form, 
we consider the approach to the bulk behavior. For suffi- 
ciently large n, one can take the continuum limit, and the 
Hessian becomes H = (e — e c (oo))9(x) — (d 2 l3/a)d 2 9/dx 2 , 
where e retains its meaning as the applied strain, d is 
the distance between adjacent layers and x is the posi- 
tion within the bulk region along the c-axis. Here, the 
layer-dependent rotation angle {—l) l 9i has now become 
the smooth function 9{x). Finally, in passing to the 
continuum limit, there is a boundary condition associ- 
ated with the interfaces, which for our simple model is 
(e c (oo) — e)9(x) = (d(3/a)d9/dx, for x = 0, nd. The criti- 
cal point where this continuum Hessian ceases to be pos- 
itive definite corresponds to the lowest eigenvalue of the 
standard one-dimensional particle-in-a-box problem of 
quantum mechanics with a modified boundary condition. 
The resulting critical strain is e c = e c (oo) — (f3q 2 )/(an 2 ), 
where q satisfies the transcendental equation tan(q>/2) = 
an(e c (oo) — e) /(J3q), whose solution approaches a con- 
stant, q — > 7r, for large n. The continuum limit thus 
reproduces exactly the 1/n 2 approach to the bulk value 
observed in the ab initio data. The slope 7 observed in 
the ab initio data can now be identified as a combination 
of parameters from the Ginzburg-Landau Hamiltonian, 
namely 7 = {j3/a)iT 2 . 

Using the observed slope in the ab initio data as well 
as data from the bulk and n = 2 phases, we now deter- 
mine completely the parameters in the Ginzburg-Landau 
model. In the bulk phase, the last term of Eq. [5] vanishes 
due to the AFD order, and the analysis becomes identi- 
cal to that for n = 1 phase, in which A, a and e c (oo) can 
all be determined from the linear fits to the 9 2 and V AE 
versus e data. When n = 2, by symmetry, the Hamilto- 
nian again assumes a form of Eq. [TJ but with e instead 
of e c (oo), so that we directly read off e = e c (n = 2). The 
final parameter remaining is /?, which we determine from 
the above result for 7 in terms of /3 and a, taking 7 from 
the slope of the best-fit line to the ab initio data (solid 
line, lower panel of Figure [3]). 

The open circles in the lower panel of Figure [3] rep- 
resent the critical strains from eigenvalues of the above 
Hessian for inhomogeneous systems, with parameters ex- 
tracted from the ab initio results as above. The clear 
linear behavior confirms our analytic prediction that the 
critical strain in the model approaches the bulk value as 
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FIG. 4: (Color online) Out-of-plane lattice constant versus 
the biaxial strain e: ab initio results (+), linear fits to data 
on either side of the transition (dashed lines). 



1/n 2 . The upper panel of Figure [3] presents an even more 
compelling case for the model free energy. It compares 
ab initio and model-Hamiltonian results for a quantity 
not included in the parameter fitting at all, the layer- 
dependent rotation angles 9i for the n = 5 RP-n phase 
at e = 0% strain. The agreement is impressive consid- 
ering the simplicity of the model and the fact that the 
model is not fit to these quantities. 

Experimental signatures — In terms of an experimen- 
tally observable signature of this transition, we find that, 
associated with the rotational ordering, there is a signif- 
icant expansion of the out-of-plane lattice constant c, by 
an amount which should be readily detectable with x-ray 
techniques. Figure 2] shows the lattice constant c of the 
n = 2 phase as a function of e through the AFD transition 
point. The transition is clearly evident where the slope 
of the curve changes substantially between e = — 1% and 
—2% strain. The lattice constant c for the AFD phase is 
significantly larger (> 0.6% for strains e < —2.6%, about 
1% strain beyond the transition) than the extrapolated 
value from the non-rotated phase. A perhaps more prac- 
tical approach is to grow various members of the RP- 



n series on substrates of fixed epitaxial strain e and to 
observe the dependence c(n) of the out-of-plane lattice 
constant as a function of n, looking again for a kink at 
the transition point. We stress that the numerical values 
given here apply specifically to zero temperature. 

Conclusion — This work presents a detailed ab initio 
study of the effects of strain on the rotational instability 
in the Ruddlesden-Popper series in strontium titanate for 
n = 1, • • • ,5. We find a similar second-order structural 
phase transition to what is seen in bulk strontium ti- 
tanate, but with a critical point displaced by ~ 1 /n 2 from 
the bulk value. The key microscopic mechanism for the 
rt-dependence is the interplay of the local AFD coupling 
between neighboring octahedra layers with a depression 
of the rotational instability for octahedra immediately 
neighboring the SrO stacking faults. It appears that the 
depression is associated with the extra freedom which 
the neighboring octahedra have to distort by formation 
of local ionic electric dipole moments and by rumpling 
of the extra SrO layers. A simple Ginzburg-Landau free- 
energy expression generalized for inhomogeneous systems 
by inclusion of nearest-neighbor inter-octahedral inter- 
actions captures the essential features of the ab initio 
results for the n-dependence of the transition. Despite 
its simplicity, this free-energy expression gives reason- 
able predictions even for quantities to which it is not fit, 
such as the distribution of rotation angles in the bulk 
regions of the RP phase. This success suggests that in- 
clusion of inter-octahedral interactions into more general 
bulk Hamiltonians, particularly those with ferroelectric 
degrees of freedom, will be a fruitful direction to pur- 
sue in future studies of nano-structured and superlattice 
perovskite materials. 
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